The BES /o(1810): a new glueball candidate 
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We analyze the /o(1810) state recently observed by the BES collaboration via radiative J/^ decay 
to a resonant (f)uj spectrum and confront it with DM2 data and glueball theory. The DM2 group 
only measured uju decays and reported a pseudoscalar but no scalar resonance in this mass region. 
A rescattering mechanism from the open flavored KK decay channel is considered to explain why 
the resonance is only seen in the flavor asymmetric uj(j) branch along with a discussion of positive C 
parity charmonia decays to strengthen the case for preferred open flavor glueball decays. We also 
calculate the total glueball decay width to be roughly 100 MeV, in agreement with the narrow, newly 
found fo, and smaller than the expected estimate of 200-400 MeV. We conclude that this discovered 
scalar hadron is a solid glueball candidate and deserves further experimental investigation, especially 
in the KK channel. Finally we comment on other, but less likely, possible assignments for this state. 



I. INTRODUCTION 

Significant in the recent wave of particle discovery, the 
BES collaboration has just reported [1] a scalar hadron 
with mass about 1812 MeV and width of 105(20) MeV. 
This /o(1810) state appeared as a 95 event enhancement 
in the Locf) spectrum accompanied by a radiative photon 
from 5.8 x 10^ J/4> decays. This paper considers several 
interpretations for this state and focuses on the most 
exciting assignment, the long-sought scalar glueball. 

Glueballs have been predicted in lattice gauge calcu- 
lations [2] and many body theory [3, 4] with both ap- 
proaches agreeing the ground state has quantum num- 
bers 0++ and mass in the range 1700 to 1800 MeV. 
Scalar hadrons between 1 and 2 GeV, predominantly the 
/o(1370), /o(1500) and /o(1710), have been scrutinized 
for glueball wavefunction components in numerous stud- 
ies [5, 6]; however, firmly identifying gluonic degrees of 
freedom remains elusive [7]. In this paper we dispel sev- 
eral theoretical conjectures about the scalar glueball and 
show that the discovered BES state is a good glueball 
candidate meriting more careful study. Using a QCD- 
based model we calculate that the full glueball decay 
width is about 100 MeV which is consistent with the 
measured 105 ± 20 MeV width for the /o(1810). We also 
show that the commonly assumed flavor blind glueball 
decay treatment entails large corrections yielding a mea- 
surable uj(j) branch but a suppressed lolo channel, again 
consistent with new data. Finally we demonstrate how 
the rescattering mechanism, /o(1810) KK — + (j)uj, fa- 
cilitates observing (puj cleanly above the tail of the pre- 
dominant /o(1710) KK spectrum, given sufficient pre- 
cision. 



II. PHENOMENOLOGICAL CONSIDERATIONS 
A. Contrasting BES with DM2 and Mark III data 

About two decades ago the DM2 collaboration ob- 
served at Orsay 8.6 x 10^ J/t/i decays and studied the 
ujuj spectrum triggered by radiative decay (photon de- 
cays have a 1% branch due to aEKi/cts suppression). 
They reported [8] a branching ratio, B{J /'ip 'yujuj) = 
1.4(0.2)(0.4) X 10-3, much larger than the BES -/Lu<j) ra- 
tio, B{J/i; -fLu<j)) = 2.61(0.27)(0.65) x lO"'^, and ob- 
served a strong pseudoscalar ?7(1760) enhancement but 
concluded there was no relevant structure near or above 
1800 MeV. Of course the ujuj decay must correlate with 
about 3 pp decays and indeed Ref. [9] observed a sim- 
ilar pseudoscalar signal in pp about 100 MeV wide, but 
threshold effects make it difficult to compare these two 
experiments. Related, there should also be a correla- 
tion with the ujip channel [10, 11] (note DM2 was not 
designed to detect uj(p) however this can only be ob- 
served for resonances above the uj(t) 1802 MeV threshold 
(it would also be kinematically suppressed unless signifi- 
cantly above threshold). 

Most recently, the BES collaboration also reported 
[12] an analysis of J/ip — > jujuj from the same 5.8 
xlO'^J/tjj BESII detector data set. They clearly confirm 
the 77(1760) (mass 1744 MeV, width 244 MeV) and also 
evidence for a 0^"*" structure which could correspond to 
the /o(1710) and/or /o(1810). However, it is difRcuh to 
determine the mass and width of this scalar hadron due 
to the dominant contributions from the 77(1760). 

We submit the absence of the /o(1810) scalar reso- 
nance, independent of its quark-glue structure, in the 
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FIG. 1: Ratio of uuj to uo(f> for a Breit-Wigner folided with 
phase space. 



DM2 an(i BESII lolo spectra iiKiicates the /o(1810) (decay 
is not flavor blind. This is because the BESll LO(j) signal, 
combined with flavor independent decay (omitting phase 
space effects), predicts an ujlj signature that would have 
also been observed by BESII and DM2. To quantify, we 
compute the ratio, i?, of lolj to LO(j) phase space factors, 
P{m), folded with a Breit-Wigner (see Fig. 1), 
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Multiplying R by the number of BES /o(1810) observed 
events (95) and the ratio of DM2 to BES J/^/- total de- 
cays (8.6 X 10^/5.8 X 10^) yields 58 /o lolj events DM2 
would have reported, assuming equal reconstruction effi- 
ciencies and a /o(1810) flavor blind decay. However, as 
detailed in Fig. 2, DM2 only reported 4 lolo events in 
the 0++ channel, which undermines the /o(1810) flavor 
blind decay assumption. Predictions for other possible 
/o mass assignments are summarized in Table I. Also, 
rescaling the BES data sample size for appropriate com- 
parison (see dashed line in Fig. 2), yields at best only a 
few (less than 4) uoff) events that the DM2 collaboration 
would be expected to observe. It would appear that the 
number of events in the DM2 J/?/; sample is insufficient 
to determine whether a resonance is or is not present in 
their Locj) spectrum even if they were looking for it. 

The DM2 data is in overall agreement with the Mark 
III data at SLAC (see [13] and references therein). The 
Mark III KK and tttt spectra featured a prominent 
/o(1710) but no /o(1810) state. Somewhat perplexing. 



FIG. 2: Dashed line: BES uocj) events rescaled to the DM2 
sample size. Solid line: actual Q^^ujlo DM2 measurement. 
A DM2 confirming ujcp decay signal is not be possible due to 
their small J /tj} sample. 



BES reported a low statistics study of the K*K* spec- 
trum in radiative J/ip decays [14] with the 0+"'" chan- 
nel not significantly populated. We also note that the 
more recent 294 events cleanly isolated by BES, 

although only part of the total produced, represent an 
extremely small branching fraction compared to, for ex- 
ample, B{J/ij -iK*K*) = 4(1) X 10^3. 

The /o(1810) decay profile is perplexing. While sup- 
pressed rates to the K*K* channel can be understood 
(note conservation of forbids decay to K*K and piv) 
due to limited phase space (see next section), what scalar 
hadron would decay leaving a clear signal in uxj) but ap- 
parently none in either of KK^ tttt or, most significantly, 
pp and LOLol In the next section we explain the sup- 
pressed pp and LOLO decays by arguing that open flavor 
(strangeness) glueball decays are favored and that KK 
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R Predicted DM2 Events 


0.105 
1.812 0.085 
0.125 


3.9 

4.1 56 
3.8 


1.793 0.103 
1.838 0.105 


5.0 71 
3.0 43 



TABLE I: iouj events DM2 would have observed in the 0++ 
channel assuming flavor blind decay. Only 4 events were 
recorded by DM2 indicating a preference for open flavor de- 
cay. Units for F and Mfg are GeV. 
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FIG. 3: Depiction of a C = + charmonium decay via glu- 
onic and quark intermediate time steps (other orderings are 
possible) . 



rescattering plays an important role in the uxj) final state. 



Radiative charmonium decay and glueball 
formation 



Consider the radiative J /ip decay to a C = +1 charmo- 
nium (on or off-shell) that subsequently decays. Having 
positive C parity favors decay via intermediate two gluon 
states and the resulting spectra shoTild therefore display 
resonances corresponding to the glueball masses. A sim- 
ple diagrammatic analysis (see Fig. 3) reveals that open 
(explicit) flavor decays, that we call "fall apart", domi- 
nate over closed (hidden) flavor decays that require color 
exchange. Here the time axis is horizontal and the (on or 
off-shell) decay sequence is: a charmed hybrid, glueball, 
light hybrid and finally a tetraquark system. This yields 
a preference for open flavor mesons (e.g. pseudoscalar 
or vector kaons) over closed (hidden) flavor uj(j) that re- 
quires final state rescattering. Table H further supports 
this point, listing established C = +1 charmonia decays 
[15] to predominantly open flavor meson states. The 
data is best interpreted by assuming a "fall-apart" decay 
mechanism with open flavor dominating over the closed 
strangeness decay that requires color exchange (rescat- 
tering). Where the (pcf) decay is unknown we enter the 
four kaon decay. Likewise where the K*K* ratio is un- 
known we listed the branching fraction to two charged 
pions and two charged kaons. This enables estimating 
upper bounds for the two-body decays. 



III. GLUEBALL DECAY WIDTHS 
A. Existing estimates 

In previous work we and others have published es- 
timates for glueball widths which we now summarize 
before presenting new computations. As detailed in 
Refs. [10, 11] the width for the decay of a scalar glue- 
ball G to two vector mesons, G — > VV , is 
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where M = 1 GeV is a fixed reference mass, Qgw 
is the GW coupling constant and k is the cm mo- 
mentum for the decay vector mesons, given hy k = 
(Mg/2)[(1 + x- x'f^ 4a;]i/2^ with x = {Mv/Mcf, 
x' = [Mv /MgY and Mq the scalar glueball mass (now 
tested against 1812 MeV). Using Vector Meson Domi- 
nance (VMD), Ref. [11] obtained gavv = 4.65, which 
gives a small partial width of 1.43 MeV for the uxj) de- 
cay reflecting the near threshold suppressed phase space. 
An independent work [16] reports a similar value for the 
coupling, gcw = 4.23. 

If we assign the larger glueball mass Mq = 1812-1-19 = 
1831 MeV, using the BES quoted upper statistical error, 
the Lucj) width increases to 7.16 MeV. For the maximum 
possible mass (including the quoted 18 MeV systemati- 
cal error) of 1849 MeV the width further increases to 15 
MeV. Table III lists predictions for other two-body de- 
cays along with parameter sensitivity. Relying on these 
results alone would suggest the glueball width is much 
broader than the /o(1810) BES candidate and its decay 
branching fraction to (jju is insigniflcant. 

Another glueball approach [17], based upon a string 
model, predicts a different, smaller total width, Tq = 140 
MeV. The decay mechanism (with inelastic rescattering 
of the kaons in the flnal state) in this constituent model 
is illustrated in Fig. 4. Their total width is much nar- 
rower than expected from VMD but more comparable 
to the width of the BES candidate. In view of this dis- 
agrcx;mc;nt bcrtween simple string and VMD estimates we 
have performed a new, first-principles calculation of the 
glueball width which is described in the next section. 



Channel 


0-+ M 0++ (x?) 1++ ixl) 2++ ixl) 


K*K* 

K+K-Tv+Tv- 
K+K*°Tr- + cc 


8.5(3) 

21(5) 5(1) 12(4) 
12(4) 32(21) 48(28) 


K+K+K-R- 


<3 

2.6(9) 1.0(4)(4) 2.4(6)(7) 
0.42(15)(12) 1.7(3)(3) 



TABLE II: Selected C = -1-1 charmonium branching fractions 
with explicit and hidden strangeness. All numbers should be 
multiplied by 10"^. 
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l^tot 


1700 4.65 

1812 

18:u 


N/P 72.1 62.8 N/P 
1.43 176 164 4.10 
7.1() 198 181 11.:! 


135 
346 
101 


1700 4.23 

1812 

1831 


N/P 59.7 52.0 N/P 
1.18 146 135 3.39 
5.92 164 153 9.38 


112 
286 
331 



TABLE 111: Partial widths (in MeV) for various values of the 
glueball mass and coupling. N/P indicates insufficient phase 
space. The final column is the sum of the other columns and 
represents a lower bound for the total width. 
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with wavefunction normalization 

dk |0(fc)|2 



(27r)3 477 

The unit normalized quark state is 



1 . 



(5) 



\ciX) = V2EY,Bic{ci)ee\n) , 



(6) 



for dressed quark creation operator B\^^ (q) and color vec- 
tor ec- Finally, the width is given by 



dr 



2Mg 



(7) 



where the four-body phase space for the final quarks is 



FIG. 4: Illustrating glueball decay with final state rescattering 
of the produced kaons to yield uxf). 



B. Ab-initio glueball width computation 

A more fundamental, QCD-based calculation for the 
total width can be obtained using many body theory [18] . 
In this relativistic, field theoretical approach an effective 
Coulomb gauge QCD Hamiltonian is approximately diag- 
onalized using the BCS and TDA many body treatments 
for the vacuum (gap equation with dressed gluons) and 
hadron states, respectively. Results are briefly described 
with further details relegated to Appendix A. 

The glueball is represented by the lightest Fock 
state consisting of two constituent BCS transverse glu- 
ons which decay to two quark pairs that subsequently 
hadronize. The decay matrix element is 



M 



dyidyHgg{-K)G{x,y)Hgg{y)\qqqq) , (3) 



where the quark-gluon Hamiltonian field interaction is 
specified in the appendix and G(x, y) is the propagator 
for intermediate scalar hybrid meson states. The scalar 

glueball state involving BCS quasi-gluon creation oper- 
ators a^", color index a, operating on the BCS vacuum 
|Q) is 



|G) = / 



dk 



(27r)3 ^ 4 



'-a'"'{k) ■ at"(-fe)|fi) , (4) 



d^i = {27r)6\MG -J^^i) ill 



dq^ 



Consult Appendix A for the remaining, technical details 
of this large-scale, multidimensional integral calculation. 
However using dimensional analysis immediately reveals 
that the total width is of order 100 MeV. Numerical pre- 
dictions are listed in Table IV. The first column is our 
reference calculation and lists the widths for a glueball 
with mass 1812 McV and a flavor independent quark- 
gluon coupling vertex. The second column is for a flavor 
dependent, and stronger, ssg vertex inspired by a Lan- 
dau gauge study in which resummed, leading Nc radia- 
tive corrections were more suppressed for light quarks. 
The dependence on flavor factors follows directly from 
Eq. (A16). The third column illustrates sensitivity to 
the glueball wavefunction. The calculation is the same 
as the first column except the TDA wavefunction is taken 
from Ref . [4] , where a slightly lighter scalar glueball mass 
of about 1725 MeV was calculated (however we maintain 
the BES 1812 MeV kinematics/phase space). Finally the 
fourth column illustrates sensitivity to the momentum 
cut-off used in the calculation and represents probability 
flux leaking to other channels. Wavefunction components 
leading to total momentum/energy above the mass of 
the decaying glueball are virtual and suppress the width. 
Eliminating them by artificially reducing the cutoff in 
the glueball wavefunction to Mq/^. increases the width 
by about a factor of 2 which is also the upper bound to 
cut-off sensitivity. 



Widths (McV) 


Flavor independent qqg vertex 


ssg ~ ^uug {ddg) 


G wf. from [4], A = 8 GeV 


G wf. from [4], A = 0.9 GeV 


Ttot 


100 


175 


50 


90 


Tlight- light 


50 


50 


25 


45 


l^light — strange 


30 


65 


15 


30 


Tstrange — strange 


15 


60 


5 


10 



TABLE IV: Total and partial widths for a scalar glueball with mass 1812 MeV. Light refers to a light u/d quark-antiquark pair 
and strange denotes a ss pair. 
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Note that in contrast to the above phcnomcnological 
models, the calculated total widths from the more funda- 
mental theory are narrower, of order 100 MeV, and con- 
sistent with the BES measurement. Indeed, our result 
also affirms an argument for narrow glueball widths pub- 
lished [19] some time ago based on the OZI rule and orig- 
inally applied to the oddball (three gluon 1~~ glueball, 
see [18]). The assertion was that charmonia decay dom- 
inantly via a glueball/oddball intermediate state, which 
in turn selects light hadron decay channels, so that the 
actual width of the glueball is about the geometric mean 
of the width of OZI- allowed and OZI-suppressed decays, 
of order a few tens of MeV. 

Concluding this subsection, our best estimate for the 
total glueball width is about 100 MeV. 

C. Width ratios including final state rescattering 

We have also evaluated final state effects and present 
further details in Appendix B. Here we qualitatively 
comment and focus on a simple rearrangement poten- 
tial between the KK and (puj channels. This flavor ex- 
change, contact potential couples different channels and 
illustrates how the uxj) signal can arise from other chan- 
nels by final state rescattering. Table V lists the ratio 
of channel 2 to 1 partial widths for different potential 
strengths. The widths are calculated using second or- 
der perturbation theory which should be reasonable as 
long as the width ratio remains below 1. Also note that 
K*K* and pp rescattering effects are suppressed by their 
large widths (50 and 150 MeV, respectively) which will 
broaden any signal and are thus not relevant to the nar- 
row BES state. More promising is the KK ujcp rescat- 
tering process which is somewhat smaller but still size- 
able. 

Other factors explain why the tttt lou) process is not 
important. Whereas the two processes 

G KK LO(p 

G —^ TTTT — > U)U> 

have similar rescattering strengths, there is a factor of 
about 2 suppression with respect to the Locp due to the 
stronger strange quark coupling, and an additional factor 
of about 4 from wavefunction overlap suppression due to 

the very different scales involved. This reduces the rela- 
tive rescattering contribution to luoj by almost an order 



1^2 


|Spin|2 |Flavor|2 \vrc\ T^/Fi 


K*K* 4>ui 
KK (t)oj 

pp — > UJUl 
TTTT — » U)U! 


25/4 2 100 MeV 0.044 
3 2 69 MeV 0.006 

25/4 9/4 106 MeV 0.020 
3 9/4 73 MeV 0.011 



TABLE V: Rearrangement potential factors and ratio of 
widths for different channels using Vss — 200 MeV. See Ap- 
pendix B for details. 



of magnitude. An /o(1810) signal in this channel would 
not be observed by DM2, and at best marginally with 
the BES statistics. 

As for the current absence of a BES /o(181G) signal 
in KK, we submit a more extensive measurement will 
observe this decay. This should include a careful exami- 
nation for any enhancement in the tail of the established 
/o(1710) ^ KK decay. Related, and as first pointed out 
in Refs. [10, 11], the ojcj) — !■ KKSir is a distinctive, novel 
glueball signature easily detected. It may be that the 
other decay channels were more difficult to observe due 
to pion background effects (e.g. pp — > An and wu) — > Gn 
and even rjr], rjr]' multiple tt). 

IV. ALTERNATIVE /o(1810) SCENARIOS 
A. Threshold cusp 

We first examine the possibility that a threshold cusp 

[20] explains the structure in the ujcf) spectrum. This kine- 
matical enhancement occurs when a two-body system in- 
elastically couples strongly to another open channel near 
threshold. Even at 1.8 GeV this condition is possible, 
and this would produce a low momentum scattering am- 
plitude having form A + B/k with k the cjcp center of 
mass momentum. However, multiplying the BES data 
by kinematical factors appropriate to each bin yields a 
resonance that is well separated from threshold which 
seems to rule out this option. 

Furthermore, the DM2 uw and pp data should have 
a similar cusp but there arc none. Rather, these data 
reveal a prominent peak, the 7/(1760), 200 MeV above 
threshold, and clearly monotonically fall towards lower 
energies. 

B. Conventional or hybrid meson 

Even though there have been many scalar meson stud- 
ies, their structure is still not completely understood. 
In the absence of mixing (claimed to be significant in 
most analyses) , quark model /o states have two isoscalar 
flavor combinations, ss and nn = {uu + dd)j\f2, with 
"^^^^Lj = ^Pq. Their ground states are shghtly above 
1 GeV [21], and for this argument we use 1.1 and 1.4 
GeV for the light and strange quark combinations, re- 
spectively. Adding 500-600 MeV for the required radial 
excitation (e.g. (/<(1020) and its radial excitation (/)(1680)) 
yields 1.6-1.7 GeV for the light, and 2 GeV for the strange 
combination. The light quark combination is marginally 
too low while the ss radial excitation state is too high to 
explain a resonance at 1.8 GeV. Moreover, one would ex- 
pect the latter to have a sizeable KK branching fraction, 
but this is not visible in the Mark III data [13] where the 
/o(1710) is dominant. Mixing the /o(1710) with a nn ra- 
dial excitation may perhaps explain the BES peak. If so 
its decay to /o(1370)7r7r might be visible, however there 
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is no simple mechanism explaining why this state should 
appear in J/-0 decays. Although this assignment cannot 
be rejected, these arguments make us suspect. 

We next examine hybrid mesons. In many body the- 
ory [22], hybrid mesons with the minimal Fock space as- 
signment qqg in an s-wavc yield a triplet (0, 1, 2)++ and 
a pseudovector 1^ . However, for typical values of the 
string tension, ^/a = 367 MeV, their masses are near but 
above 2 GeV. Similar, though a bit lighter, results are 
obtained in the flux tube model and lattice gauge the- 
ory, so one cannot discard a hybrid state as low as 1.8 
GcV. However qualitatively comparing hybrid and glue- 
ball total width calculations, the hybrid width will be 
much broader since there is only one gluon-quark ver- 
tex interaction, instead of two, yielding one less factor of 
ttg suppression. Therefore a broader state than the BES 
result is expected. 



C. The /o(1710) tail 

The mass and width of the /o(1710), another glueball 
candidate, are poorly determined and the PDG values 
are M = 1714(5) MeV, T = 140(10) MeV. However, con- 
sistent with recent BES data, these values could bo as 
high as M = 1740 MeV, T = 166 MeV. In this case, 
an overlap with the trailing edge of the /o(1710) Breit- 
Wigner distribution could produce the observed iu'<f> sig- 
nal. However this possibility appears unlikely considering 
the near-threshold behavior of the BES u}(j) spectrum. As 
in the cusp hypothesis, the current data seems to indicate 
that the resonance is separated from the threshold and 
therefore cannot stem from the /o(1710). Higher sample 
count studies would be very useful. 



D. Four quark states 



V. SUMMARY AND CONCLUSIONS 

In this work we have examined and compared indepen- 
dent J/tp decay data sets in the 1800 MeV mass region. 
Based on the data reported by the BES collaboration, 
we believe that the newly found /o(1810) is a promising 
glueball candidate or a state with a large glueball compo- 
nent. Significantly, its mass and quantum numbers arc in 
agreement with previous theoretical expectations and its 
somewhat surprising narrow width of order 100 MeV is 
consistent with new fundamental calculations. We have 
addressed the perplexing issue of its selective decay to 
the W(/> channel and discussed why it was not observed in 
both the BESII and the smaller sample DM2 measure- 
ments in the lulo channel. Also the tttt and KK channels 
have been examined by MARK HI and only the broad 
/o(1710) structure is apparent. However, the binning of 
these data is somewhat coarse and further structure can- 
not be ruled out. The DM2 data [26] shows a falling 
slope probably due to the /o(1710) tail. Higher precision 
studies are clearly needed. 

We also noted that the rescattering process, J/tl; 
jG -yKK 7^^, may be producing the BES sig- 
nal. In view of the important theoretical implications of 
a glueball state, we submit that a precision study of the 
KK spectrum is crucial and may reveal a significant en- 
hancement around 1812 MeV. Should this test fail, we 
would be facing a new theoretical puzzle. 
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Tetraquark systems, another actively investigated 
area, also appear naturally as an intermediate step in 
a J/ip decay chain. However, as in the hybrid case, one 
expects a four quark state to decay with a broad width 
generating a background, not a sharp signal, for radiative 
J/ip decay. Related, a realistic tetraquark width predic- 
tion also requires including KK rescattering effects since 
the uxp attraction is not as strong as in KK where an- 
nihilation diagrams provide attractive forces. The quark 
rearrangement coupling between the K*K* and ujcj) chan- 
nels also provides attraction (see Appendix B). This fol- 
lows from the Resonating Group Method (RGM) [23-25] 
which predicts an increased attraction between mesons 
when each has a quark and anti-quark of the same fla- 
vor. Hence if the BES state is not a glueball, the RGM 
coupled channels will play an important role in elucidat- 
ing its structure and applications of our model to this 
system are in progress. 
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APPENDIX A: COMPUTATION OF THE 
INCLUSIVE DECAY WIDTH 

Here we present the many body effective QCD Hamil- 
tonian calculation for the glueball decay to four quarks 
that subsequently hadronize. In the Coulomb gauge the 
effective quark fSamiltonian contains an instantaneous in- 
teraction, mediated by the infrared enhanced Coulomb 
potential, and a transverse gluon exchange interaction 
that is infrared suppressed via the generation of a mass 
gap [27]. First the instantaneous interaction is diago- 
nalized to obtain the glueball bound state wavefunction. 
Then the triple quark-gluon coupling interaction 



Hqg=9 j 



dx^T^ot * • A" 



(Al) 



is treated perturbatively to compute the decay ampli- 
tude. Omitting the momentum conserving delta func- 
tions arising from the commutators, the integrand, / of 
the matrix element in Eq. (3) reduces to 



I{G qmqaqi) = CT 

-m 



2u{k) 



(A2) 



V2M2^i2£22£32S4^G5c,c.C3C4(qiq2q3q4) • 



47r 

For total (inclusive) decay the color tensor is 

Ob a ■ 



78 



C1C2 C3C4 



(A3) 



which when squared and summed over the quark color 
indices Cj yields the squared color factor CF^ = 1/4. 
The above result is for only one specific flavor and below 
we include the modification for application to three light 
flavors (u, d, s). The \/2 factor is a result of gluon ex- 
change symmetry and the glueball normalization in Eq. 
(4). The gluon self-energy, uj{k), follows from the in- 
termediate gluon propagators in Fig. 3 and is the solu- 
tion of a mass gap equation that is well approximated 

by y^m^~i-"fc^ (used here). The tensor S depends on the 

spinors in the Fourier expansion of the quark field ^I* and 
the Dirac a matrices coupled to the gluon spins. These 
spinors are usually expressed in terms of a BCS angle, 
whose relation to the running mass (from the quark gap 
equation) is sin(/)(g) — Sq — m{q) / m{q)^ + (here we 
fix m(g) = m). Squaring the matrix element and sum- 
ming over spins in the final state, we find, in terms of 
unit momentum vectors 

^ |5|2 = (l + SiS2)(l + S3S4) (A4) 

C1C2C3C4 

+ (1 + siS2)c3C4k-q3k-q4 + (1 + .S3.S4)ciC2k-qik-q2 



-C1C2C3C4 



qiq3q2q4 + q2q3qiq4 qiq2q3q4 

-k q2k q4qi q3 k qik q4q2 qs + k qik q2q3 •q4 
-hk q3k q4qi q2 - k q2k q3qi •q4 - k qik q3q2 •q4 

-(-2kqikq2kq3kq4 . 



Since all possible relative angles in the final state appear, 
J is only invariant under rigid rotations of the five vectors. 
With this, the squared decay matrix element summed 
over the final state color, spin and flavor indices reads 



R 



FF^i? 



CF 



4a;(fc)2 



2M2E12E22E32E4 



{A-Kas{k)Y\GSr2SiV ■ (A5) 



Symmetry considerations apply if both emitted quark- 
antiquark pairs are indistinguishable and we absorb this 
into the flavor factor FF^ below. 

Let us now examine the phase space integrals. Momen- 
tum conservation for a decaying glueball at rest requires 



qi + q2 + q3 + q4 = , 

yielding for the gluon momentum 

k = qi + q2 = -q3 - q4 



(A6) 



(A7) 



There are nine integration variables in Eq. (8). We 
can arbitrarily fix k along the third axis and obtain a 
47r factor for global rotations. Further we can integrate 
one azimuthal angle (say 0i) around this fixed axis, an 
operation preserving all relative angles. The modulus 
k = |k| remains an independent variable. The others 
can be chosen as gi = |qi|, coaOq-^k (which automatically 
fixes qi and q2), and the three spherical coordinates of 
q3) 93 = Iqalj cos^qgfc and ^3. Only five of these six are 
independent since we have not utilized the energy con- 
servation relation. This imposes cumbersome restrictions 
on the angular variables, so it is convenient and custom- 
ary to introduce an auxiliary variable, £, representing the 
energy of the second pair, by means of 

5{MG-Y^Ei) = j 5{MG-£-Ei-E2)S{£-E3-E4)d£ . 

The two 6 functions can be used to integrate over the two 
polar angles cos Oq-^k and cos Oq^k, leaving only the £ inte- 
gration with integration limits fixed by the requirement 
that the cosine values remain in the interval (—1,1). This 
is easily implemented in our 5 dimensional Monte Carlo 
computation by rejecting points exceeding this bound. 
The resulting polar cosines are 

COS.,.. = ml + ,l + k^-iMa-£ + E.r 



cos Oq^k 



2kqi 

{£-E3r-{ml + ql + k^) 
2kq3 



(A9) 



Note that the change of variable from energy to angle in 
each of the 5 functions adds an extra factor 



6{Eo-E) = 6{Eo - y/m'^ + k'^ + q^ + 2kqcose) 

Eq 



kq 



6 (cos 60 — cos 9) 



(AlO) 
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There are then four remaining integration variables k, 
Qij 93 7 03 J for a total of five integrals that are performed 
numerically. A representation for the coupling as in the 
infrared is needed and we use [28] 



- 



91og((fc2 + M2)/A2) 



(All) 



with A ~ 0.2-0.21 GeV and Mq ~ 1-1.1 GeV. The final 
ingredient is the propagator G{x, y) for the intermediate 
hybrid meson cut in Fig. 3, necessary for a second order 
calculation. Its exact energy eigenfunction expansion is 



h= hybrid 



1 



Mc-Eh- ie 



-{h\ . 



(A12) 



The spectrum of hybrid mesons has been studied with 
this many body method in [22] where, for string tension 
\/a- = 367 MeV, the ground state scalar hybrid meson has 
mass 2100 MeV. Excitations thereof appear with spac- 
ings similar to those in ordinary meson quark models. 
We use G = 1/(Mg - Mu) -1/(300 MeV) and dressed 
quark masses m„ = rud = 100 MeV, TOs — 200 MeV, 
consistent with prior work using the same approach and 
parameters. These values are typical of the masses calcu- 
lated in Ref. [29] , but somewhat low compared to quark 
model phenomenology since in field theory approaches a 
sizeable fraction of the hadron mass originates from the 
self-energy contribution in the bound state problem and 
not in the mass gap equation. These values also yield 
a realistic conventional hadron spectrum. Similarly we 



J 



have m,g = oj{0) = 650 MeV, where uj{k) is the solution 
of the gluon gap equation of pure gluodynamics [3] . 

Finally let us examine the flavor factors. For an inclu- 
sive decay we can separate the sum over the final states 



(A13) 



flavor 



diff 



where 



FF oc {n\{uu + dd+ ss) ^ \q'q'} (A14) 



Q Q 



(g Y I {uu + dd+ ss) \ qqqq) 



We have 



FF' 



16 • 3 = 24 



(A15) 



where the first term accounts for the case where the out- 
going quark pairs have different flavors, and the second 
for the case where the outgoing flavors are the same. The 
3 in each term reflects the number of distinct choices 
for three flavors (u,d,s), and the 1/2^ corrects for over- 
counting in the sum over final states with two pairs of 
identical particles. If we separate by flavor channel the 
corresponding factors arc 12 for light- light, 8 for light- 
strange and 4 for strange-strange. With this we obtain 
the complete expression for the glueball width (see Ta- 
ble IV for numerical results) 



r : 

E2 E4 



(27r)(47r) 



Mg /-Stt 

2Trd£ 



nMG/2 /.Mg/2 

. ,w . I k dk 

(27r)y Jo Jo Jo JO 

1 2 \,p{k)\^ 



PZTT p-MgI'^ I'lVlG/'^ /•Mg/'^ 

j d4)s / k'dk I qldqi / q^dqs 
Jo Jo Jo Jo 

|5|2(47ra,(A;))2e(cos2 9g,k - l)e(cos2 Og.k - 1) ■ 



■ kqs kqi {Mg - Mhf 



CF^FF^ 



(A16) 



4<^fc 47r 



APPENDIX B: RESONATING GROUP METHOD 
AND DECAY CHANNEL RECOUPLING 



1. Coupled channels 



In this appendix we theoretically treat the sequential 

decay of a glueball G to a meson pair followed by rcar- 



We approximately solve the equation of motion, _ff $ 



rangement. To be speciflc, we assume that G first decays E<^, using the resonating group/coupled channels formal- 



to K*K* and then to (fxjo as depicted in Fig. 4. 



ism [23-25] for this three channel problem 



I Hg-E 

\ 



Vs 



sb 

V* 

re 





ie V, 
I 



re 




= 



(Bl) 



where Vsb is the string breaking decay coupling between arrangement potential coupling the latter to the (jxjj chan- 
the glueball and the open flavor channel, and Vre is the re- 
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ncl. Wc now extract the glucball width F from the imag- 
inary part of the resuhing ghieball energy/mass. Since 
Hq — Mg = 1812 McV and from our computation of 
the total glueball width, it follows that Vsb is at most 
of order 100 MeV. We can therefore diagonalize using 
perturbation theory to the leading order in Vab, 



The only model dependent quantities arc the poten- 
tial strengths Vre and Vsi, but only one enters the ra- 
tio of the two partial decay widths. Because the vector 
mesons have the same oscillator parameter and the re- 
duced masses Hc are similar, while the thresholds Mg 
differ for K*K* and (buj, the ratio reduces to 



Ha-E-V, 



sb 



1 



1 



V* 



Hi^.fi.-E-ie H^^-E-ie 
1 



Hk'K' - E-ie 



sb 



$G = 



H,ff{E)--T{E)-E]^G = 0. (B2) 



Also using perturbation theory for Vre, we identify the 
imaginary part of the first potential term, to order Vj,, as 
the partial decay width to K*K*, while the contribution 
(for an open (j)oj channel) to the imaginary part from the 
second term, to order Ygf^V^e, yields the partial decay 
width for 0a;. The same analysis can be used for other 
sequential decays, e.g. G — > tttt — > wu. 

In evaluating the string breaking and rearrangement 
potentials, we truncate the sum over intermediate hadron 
states to the ground state mesons in each channel c and 
employ harmonic oscillator wavefunctions, ^q'^, having 
oscillator parameter Uc- This yields the following sepa- 
rable potentials for string breaking, 



and rearrangement, 

Vre = Vre\(po''){^0°\ ' 

and resulting partial widths 

r^.^. = 2\vsb\^Im[gK,K,iE)] , 

^4,0, = 2\Vsb\'^\Vre\'^Re[gK-'if{E)]>^ 

Im [g^u^{E)] Re [qk'k- {E)] , 



(B3) 
(B4) 



9ciE) = (0?' 



1 



Mr 



E 



I0«^). (B5) 



Here Me and /Xc are the threshold energy and reduced 
mass for channel c. Because the glueball mass is near 
the channel thresholds, the real and imaginary parts of 
the channel Green functions, gc{E) = Uc + ibc, can be 
well approximated by 



Jo 



-q'^dq 



be ~ 4ar^^le^/2^^fle^/E - Me 
and E = Mg- The partial decay widths are then 



K'K' 



2\Vsb\'^\Vre\'^a'l.. j^.b^u: 



(B6) 
(B7) 

(B8) 



/ . 2 \2 / Mg-m^^ 

= [Vre4.H4u,a^J ' 



Mg - Mk^k, 



■ (B9) 



The parameters used in our calculation are listed in Ta- 
ble VI. Notice that we did not compute the complete 
geometric series in Eq. (B2). If the ratio, Eq. (B9), 
is large, we need to sum the full geometric series in Eq. 
(B2). However, the remaining terms of the series con- 
tribute the same factor for the decay to K*K*, or to (f)co, 
and therefore the ratio is correct to all orders in Vre- 



2. String breaking 

Because the ratio in Eq. (B9) simplifies, we only list 
the flavors directly produced with string breaking. Notice 
that the same flavors are produced with a direct decay 
of the constituent gluons. Let us suppose that there are 
two string breakings producing two mesons. We assume 
that each string breaking creates a quark-antiquark pair 
in an approximately symmetric way, yielding an SU(3) 
flavor singlet 



dd 



(BIO) 



where we suppress spin and color notation. In the string 
breaking picture, the quarks will be separated, each to 
one of the two produced mesons. This is also necessary 
to ensure that each of the produced mesons is a color sin- 
glet (equivalent to having a quark exchange between the 
two flavor singlet sources and the two produced mesons). 
Exchanging the first and third quarks (permutation op- 
erator P^^) yields 



3l3| 



{uu + dd-\ 
f dddd - 
+USSU + suus 



= \uuuu 



ss)(uu -\- dd-\- ss)) 
- uddu + duud 
f dssd + sdds + ssss) 



(Bll) 



uu + dduu + dd uu — dduu — dd 

+uddu + duud+ ussu + suus + dssd + sdds + ssss) 

with a similar result for the exchange of the second and 
fourth antiquarks {P^'^)- 

Specializing to vector- vector production, nine different 
vector pairs are produced and Eq. (Bll) becomes 

(B12) 



I WW + p°p° + p'^p~ + p~p'^ 



+K*+K* 



K*-K* 



but no u><p from string breaking. This explains why the 
off-diagonal Hamiltonian matrix element coupling this 
channel in Eq. (Bl) is zero and that uxp can only be 
produced by rearrangement. 
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Rearrangement 



so that 



In evaluating T^^^ for the partial width ratio, Eq. (B9), 
the normalized \K*K*) and \u}^) states are needed. The 
oj<p component is 



dd 



V2 



ss) , 



(B13) 



while K*K* is given by 

X*'^K*~ + K*~ K*~^ + K*'^K*^ -\- j{*oj^*o 



\K*K*) 



V2V4 

uddu + duud + ussu + suus 



(B14) 



where the wavefunction normalization includes the 
meson-meson exchange, 



{K*K*\l + P^^P'^'^\K*K*) = 1 . 



(B15) 



Then the flavor rearrangement matrix element involving 
P^^ quark-quark exchange is 



{K*+K*-\P''\<Pu) 



(B16) 



Notice that P^^ antiquark-antiquark exchange produces 
exactly the same result. This also applies to the color and 
spin X space rearrangement overlaps, so we only compute 
the P^^ overlaps and include an additional factor of 2 to 
account for antiquark exchange. 
The color rearrangement overlap is 



(1 1|P^^|1 1) = - 



(B17) 



For the space x spin matrix clement, using the graphi- 
cal rules [30, 31], a separable potential with strength v 
emerges 



{cj>Zci>z\P^^\4^Z'i>o) = v\cj>-^^)m 

Hence, only the hyperfinc interaction, 
A,; • A, 



Vss 



-16/3 



Si ■ Sj 



(B18) 



(B19) 



contributes to the space x spin rearrangement and the 
potential v is proportional to Vss — 200 MeV obtained 
from Ref. [32]. To determine this constant we consider 
the specific overlap where both K*K* and (fxjo have a 
total spin 0. Coupling two vector mesons yields 



jOO) 



111)11 - 1) - |i())|io) + |i - 
TTii + iiTT nn + TiiT + itn + iTiT 



^/3 



2^3 



(B20) 



(00|pi3|00) = . (B21) 

Only the intra-cluster contributions V13, V14, V23, V24 
of the hyperfine potential need to be considered [33, 34] 
since the inter-cluster ones are already included in the 
meson mass calculation. Adding all intracluster contri- 
butions, the total space x spin overlap contribution for 
this case is 



-7:Vss\<t>^){^^\ . 



(B22) 



Then the resulting color x flavor x space x spin over- 
lap contribution for this case is 



(B23) 



yielding the strength and sign for the rearrangement po- 
tential 



1, 



':Vss ■ 



(B24) 



The value of a is fixed by the rms radius, < >, of the 
corresponding meson. For a Gaussian wavefunction 



< > = 



dr r^e-'-'/"' 



/o 
3 , 



(B25) 



For the pion < >^/^~ 0.5 fm, but this is anomalously 
large due to the light mass (in chiral perturbation theory 
it is divergent in the chiral limit). We use 0.4 fm for 
kaons and 0.3 fm for all vector mesons. Converting to 
MeV~^, Table VI lists the oscillator parameters used for 
each channel. 



c — > 


4>uj ujui K*K* pp KK 'K-n 


lO^ae(MeV-i) 
10^&c(MeV-i) 
a-\MeV) 
Mc(McV) 
Mc(MeV) 


2.74 2.41 2.75 2.38 2.16 1.18 
0.57 2.35 0.96 2.40 3.51 1.20 
804 804 804 804 603 483 

443 391 446 385 248 69 
1802 1564 1784 1540 994 278 



TABLE VI: Rescattering parameters. 
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